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ABSTRACT 

A large international effort is under way to assess the presence of a shadow in the radio emis¬ 
sion from the compact source at the centre of our Galaxy, Sagittarius A* (Sgr A*). If detected, 
this shadow would provide the first direct evidence of the existence of black holes and that 
Sgr A* is a supermassive black hole. In addition, the shape of the shadow could be used to 
learn about extreme gravity near the event horizon and to determine which theory of gravity 
better describes the observations. The mathematical description of the shadow has so far used 
a number of simplifying assumptions that are unlikely to be met by the real observational 
data. We here provide a general formalism to describe the shadow as an arbitrary polar curve 
expressed in terms of a Legendre expansion. Our formalism does not presume any knowledge 
of the properties of the shadow, e.g., the location of its centre, and offers a number of routes to 
characterize the distortions of the curve with respect to reference circles. These distortions can 
be implemented in a coordinate independent manner by different teams analysing the same 
data. We show that the new formalism provides an accurate and robust description of noisy 
observational data, with smaller error variances when compared to previous approaches for 
the measurement of the distortion. 

Key words: black hole physics - Galaxy; centre - submillimetre: galaxies 


1 INTRODUCTION 

There is a widespread belief that the most convincing evidence for 
the existence of black holes will come from the direct observa¬ 
tions of properties related to the horizon. These could be through 
the detection of gravitational waves from the collapse to a rotating 
star (Baiotti & Rezzolla 2006), from the ringdown in a binary black 
hole merger (Berti et al. 2009), or through the direct observation 
of its ‘shadow’. In a pioneering study, Bardeen (1973) calculated 
the shape of a dark area of a Kerr black hole, that is, its ‘shadow’ 
over a bright background appearing, for instance, in the image of 
a bright star behind the black hole. The shadow is a gravitationally 
lensed image of the event horizon and depends on the closed or¬ 
bits of photons around the black hole^. Its outer boundary, which 
we will hereafter simply refer to as the shadow, corresponds to the 
apparent image of the photon capture sphere as seen by a distant 
observer. General relativity predicts, in fact, that photons circling 
the black hole slightly inside the boundary of the photon sphere 

^ Strictly speaking, also an horizon-less object such as a gravastar (Mazur 
& Mottola 2004) would lead to a shadow. However, rather exotic assump¬ 
tions on the heat capacity of the gravastar’s surface are needed to jus¬ 
tify a lack of emission from such a surface (Broderick & Narayan 2007); 
gravitational-wave emission would unambiguously signal the presence of 
an event horizon (Chirenti & Rezzolla 2007). 


will fall down into the event horizon, while photons circling just 
outside will escape to infinity. The shadow appears therefore as a 
rather sharp boundary between bright and dark regions and arises 
from a deficit of those photons that are captured by the event hori¬ 
zon. Because of this, the diameter of the shadow does not depend on 
the photons’ energy, but uniquely on the angular momentum of the 
black hole. In general relativity and in an idealized setting in which 
everything is known about the emission properties of the plasma 
near the black hole, the shadow’s diameter ranges from 4.5 for 
an extreme Kerr black hole, to fs for a Schwarzschild black 
hole, where := 2GM/(? is the Schwarzschild radius. In prac¬ 
tice, however, the size and shape of the shadow will be influenced 
by the astrophysical properties of the matter near the horizon and, 
of course, by the theory of gravity governing the black hole. 

Besides providing evidence on the existence of black holes, 
the observation of the black hole shadow and of the deformations 
resulting in the case of nonzero spin, is also expected to help de¬ 
termine many of the black hole properties (see e.g. Chandrasekhar 
1998; Falcke & Markoff 2013; Takahashi 2004; Falcke et al. 2000; 
Doeleman et al. 2009). More specifically, imaging the shadow of a 
black hole via radio observations will allow one to test the predic¬ 
tions of general relativity for the radius of the shadow and study 
astrophysical phenomena in the vicinity of black holes [see Jo- 
hannsen & Psaltis (2010) and, more recently, Zakharov (2014) for 
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a review]. In addition, it will allow one to set constraints on the 
validity of alternative theories of gravity which also predict black 
holes and corresponding shadows (see e.g. Eiroa & Sendra 2014; 
Tsukamoto et al. 2014; Falcke & Markoff 2013). 

The possible observation of a black hole shadow has recently 
received a renewed attention as the spatial resolution attainable 
by very long baseline interferometry (VLBI) radio observations is 
soon going to be below the typical angular size of the event horizon 
of candidate supermassive black holes (SMBHs), such as the one 
at the centre of the Galaxy or the one in the M87 galaxy (Falcke 
et al. 2000). These observations are the focus of international sci¬ 
entific collaborations, such as the Event Horizon Telescope (EHT)^ 
or the Black Hole Camera (BHC)^, which aim at VLBI obser¬ 
vations at 1.3 mm and 0.87 mm of Sagittarius A* (Sgr A*) and 
M87. We recall that Sgr A* is a compact radio source at the cen¬ 
tre of the Galaxy and the SMBH candidate in our galaxy. In fact, 
the orbital motion of stars near Sgr A* indicates that its mass is 
~ 4.3 X 10®Mq (Ghez et al. 2008; Genzel et al. 2010). 

Given a distance of 8 kpc from us, the angular size of the 
Schwarzschild radius of the SMBH candidate in Sgr A* is ~ 
10 /ras, so that the corresponding angular diameter of the shadow 
is of the order of ~ 50 /ras. Similarly, with an estimated mass of 
~ 6.4 X 10® Mq (Gebhardt et al. 2011) and a distance of 16 Mpc, 
the M87 galaxy represents an equally interesting SMBH candidate, 
with an angular size that is of the same order i.e., ~ 40 /ras. Al¬ 
though the resolution achievable at present is not sufficient to ob¬ 
serve an image of the shadow of either black hole, it is sufficiently 
close that it is realistic to expect that near-future observations will 
reach the required resolution. Indeed, future EHT and BHC obser¬ 
vations of Sgr A* are expected to go below the horizon scale and to 
start to provide precise information on the black hole orientation, as 
well as on the astrophysical properties of the accretion flow taking 
place on to the black hole (Chan et al. 2015; Psaltis et al. 2015). 

An extensive literature has been developed to calculate the 
shadow of the black hole in known space-times, either within gen¬ 
eral relativity (Young 1976; Perlick 2004; Abdujabbarov et al. 
2013), or within alternative theories of gravity (Bambi & Modesto 
2013; Tsukamoto et al. 2014; Amarilla et al. 2010; Amarilla & 
Eiroa 2012, 2013; Atamurotov et al. 2013a,b; Schee & Stuchlfk 
2009). In most cases, the expression of the shadow as a closed polar 
curve is not known analytically, but for the Plebahski-Demiahski 
class of space-times, the shadow has been cast in an analytic 
form (Grenzebach et al. 2014, 2015). 

Because the shadow is in general a complex polar curve on 
the celestial sky, an obvious problem that emerges is that of the 
characterization of its deformation. For example, in the case of a 
Kerr black hole, the difference in the photon capture radius be¬ 
tween corotating and counter-rotating photons creates a “dent” on 
one side of the shadow, whose magnitude depends on the rotation 
rate of the black hole. A way to measure this deformation was first 
suggested by Hioki & Maeda (2009) and then further developed 
by other authors (Bambi & Freese 2009; Bambi & Yoshida 2010; 
Bambi & Modesto 2013). In essence, in these approaches the shape 
of the shadow is approximated as a circle with radius Rg and such 
that it crosses through three points located at the poles and at the 
equator of the shadow’s boundary. The measure of the dent is then 
made in terms of the so called ‘deflection’, that is, the difference 
between the endpoints of the circle and of the shadow, with a di¬ 


® http://eventhorizontelescope.org/ 
^ http://blackholecam.org/ 


mensionless distortion parameter being given by the ratio of the 
size of the dent to the radius Rs [cf. Eq. (48)]. 

While this approach is reasonable and works well for a black 
hole such as the Kerr black hole, it is not obvious it will work 
equally well for black holes in more complex theories of gravity 
or even in arbitrary metric theories of gravity as those considered 
by Rezzolla & Zhidenko (2014). Leaving aside the fact in all these 
works the shadow is assumed to be determined with infinite preci¬ 
sion (an assumption which is obviously incompatible with a mea¬ 
sured quantity), many but not all approaches in characterizing the 
black hole shadow and its deformations suffer from at least three 
potential difficulties. They often assume a primary shape, i.e., that 
the shadow can be approximated with a circle; exceptions are the 
works of Kamruddin & Dexter (2013) and Psaltis et al. (2014). 
In the first one, a model has been proposed to describe the “cres¬ 
cent” morphology resulting from the combined effects of relativis¬ 
tic Doppler beaming and gravitational light bending (Kamruddin 
& Dexter 2013); in the second one, an edge-detection scheme and 
a pattern-matching algorithm are introduced to measure the prop¬ 
erties of the black hole shadow even if the latter has an arbi¬ 
trary shape (Psaltis et al. 2014). (ii) They assume that the observer 
knows the exact position of the centre of black hole^. (Hi) They 
are restricted to a very specific measure of the distortion and are 
unable to model arbitrary distortion; exceptions are the works of 
Johannsen & Psaltis (2010) and Johannsen (2013), where instead 
polar-averaged distortions have been proposed. 

To counter these potential difficulties, we present here a new 
general formalism that is constructed to avoid the limitations men¬ 
tioned above. In particular, we assume that the shadow has an ar¬ 
bitrary shape and expand it in terms of Legendre polynomials in a 
coordinate system with origin in the effective centre of the shadow. 
This approach gives us the advantage of not requiring the knowl¬ 
edge of the centre of the black hole and of allowing us to intro¬ 
duce a number of parameters that measure the distortions of the 
shadow. These distortions are both accurate and robust, and can 
be implemented in a coordinate independent manner by different 
teams analysing the same noisy data. 

The paper is organized as follows. In Sect. 2 we develop the 
new coordinate-independent formalism, where an arbitrary black 
hole shadow is expanded in terms of Legendre polynomials. Using 
this formalism, we introduce in Sect. 3 various distortion parame¬ 
ters of the shadow. In Sect. 4 we apply the formalism to a number 
of black hole space-times by computing the coefficients of the ex¬ 
pansion and by showing that they exhibit an exponentially rapid 
convergence. We also compare the properties of the different dis¬ 
tortion parameters and assess which definition appears to be more 
accurate and robust in general. Section 5 offers a comparison be¬ 
tween the new distortion parameters introduced here with the more 
traditional ones simulating the noisy data that are expected from 
the observations. Finally, Sect. 6 summarizes our main results the 
prospects for the use of the new formalism. 

We use a system of units in which c = G = 1, a space¬ 
like signature (—,-I-,and a spherical coordinate system 
{t, r, 6, (j)). Greek indices are taken to run from 0 to 3. 


^ The reconstruction procedures of the observational data does resolve this 
problem, which is however still present when considering and comparing 
purely theoretical representations of the shadow 


Characterization of a black hole shadow 3 



Figure 1. Schematic representation of the biack hoie shadow as a generic 
poiar curve in a coordinate system («, /3) with origin O in the ‘centre’ 
of the shadow. The iatter is transiated by a vector Ro with respect to the ar¬ 
bitrary coordinate system (o', /?') with origin O' in which the observations 
are made. 


2 GENERAL CHARACTERIZATION OF THE SHADOW 

In what follows we develop a rather general formalism to describe 
the black hole shadow that radio astronomical observations are ex¬ 
pected to construct. For all practical purposes, however, we will 
consider the problem not to consist of the determination of the 
innermost unstable circular orbits for photons near a black hole. 
Rather, we will consider the problem of characterizing in a mathe¬ 
matically sound and coordinate-independent way a closed curve in 
a fiat space, as the one in which the image will be available to us as 
distant observers. 

Assume therefore that the astronomical observations provide 
the shadow as an one-dimensional closed curve defined by the 
equation 


R' = R'{'b'), ( 1 ) 

where R! and ip' can be thought of as the radial and angular coor¬ 
dinates in a polar coordinate system with origin in O'. In practice, 
astronomical observations will not be able to provide such a sharp 
closed line and a more detailed analysis would need to take the ob¬ 
servational uncertainties (which could well be a function of ip') into 
account. We will discuss some of these uncertainties in Sect. 5, but 
for the time being we will consider the shadow as an idealized one¬ 
dimensional closed curve. A schematic example of the polar curve 
is shown in Fig. 1, where a' and f are the so-called “celestial co¬ 
ordinates” of the observer, and represent an orthogonal coordinate 
system with one of the unit vectors being along the line of sight. 

Of course, there is no reason to believe that such a coordinate 
system is particularly useful, or that in using it a nonrotating black 
hole will have a shadow given by a perfect circle. Hence, in order 
to find a better coordinate system, and, in particular, one in which 
a Schwarzschild black hole has a circular shadow, we define the 
ejfective centre of the curve extending the definition of the position 


of the centre of mass for a solid body to obtain 


Rq ~ 


Jo 


e„,(ip')R' [gR'it'(dR'/dip'f 


11/2 


dip' 


Jo"" [9R'R’{dR'/dip'f dip' 


( 2 ) 


where e^, is the radial-coordinate unit vector and where gwRi, 
are the metric functions of the polar coordinate system 
{R', Ip'). Two important remarks: first, radio observations may well 
yield, especially in the nearest future, only a portion of the shadow, 
namely, the one with the largest brightness. Yet, it is useful to con¬ 
sider here the shadow as a closed polar curve since this is the way it 
is normally discussed in purely theoretical investigations. Second, 
at least observationally, the location of the centre of the black hole 
shadow is a free parameter in the image-reconstruction procedure 
and so already part of the analysis of the observational data. At the 
same time, the definition of a centre is useful also in the absence 
of actual observational data since it can help in the comparison of 
shadows that are built analytically and hence without observational 
data. 


From the knowledge of the vector Ro, the coordinate position 
of the effective centre can be expressed explicitly in terms of the 
radial and angular coordinates as 


Ro ■= 


/ ^27r \ ^ / /’27r 

:= R'dip'j R''^ COS ip'dip' 


('2'k \ 2 

R'^ sin Ip'dip' 


1/2 


ipo ~ tan 


f " R'^ sin Ip'dip' 
R'^ cos Ip'dip' 


(3) 

(4) 


We note that if the centre of the primary coordinate system O' 
coincides with the black hole origin, then the parameter Ro exactly 
corresponds to the shift of the centre of the shadow with respect to 
the black hole position defined by Tsukamoto et al. (2014). 

Having determined the effective centre of the shadow, it is 
convenient to define a new polar coordinate system centred in it 
with coordinates {R, ip). Clearly, the new coordinate system with 
origin O is just translated by .Rq with respect to O' and, hence, the 
relation between the two coordinate systems is given by 


i? := (R' cos Ip' — Ro cos ipof 

.I 1 1/2 

+ (R'sint/i'— Ro sint/io) , (5) 

_i R'sin t/i'- Rosin t/)o 

Ip—tan —- - - - - -. (6) 

R' cos Ip' — Ro cos Ipo 


Note that we have kept the new axes a and j3 parallel to the origi¬ 
nal ones a' and pt'. This is not strictly necessary but given the ar- 
bitrarity of the orientation of both sets of axes, it provides a useful 
simplification. 

A well-defined centre of coordinates allows us now to obtain 
a robust definition of the reference areal circle as the circle hav¬ 
ing the same area as the one enclosed by the shadow. In particular, 
given the closed parametric curve R = R(V')> its area will in gen¬ 
eral be given by 



where in the second equality we have set gRR = 1, = R{ip), 
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Figure 2. Schematic representation of the local distortion between the 
polar curve representing the black hole shadow and representative cir¬ 
cles with circumference and area radii and R ^, respectively. 


while in the third equality we consider the representation of the 
curve in terms of a more generic parameter A, i.e., R = R{tp{X)). 
If the shadow is a closed curve, the integration limits Ai _2 can be 
found from the condition ’tp{X) = 0 and V’(A) = 2tt, while they 
will be restricted by the actual observational data when only a por¬ 
tion of the shadow is available. 

We can then define the areal radius R^ of the reference circle 
simply as 


R 



( 8 ) 


Similarly, and if simpler to compute, it is possible to define the 
circumferential radius Rc of the reference circle as 


• 27r ’ 

where the circumference is calculated as 


C~J {dR^ + 



(9) 


( 10 ) 


3 DISTORTION PARAMETERS 


With a well defined and unambiguous set of coordinates {R, tp), we 
can next move to the characterization of the geometrical properties 
of the shadow. To this scope we simply employ an expansion in 
terms of Legendre polynomial, i.e., we define 


;= ^qP^(cos?/^) , (13) 

£=0 


where Pi{costp) is the Legendre polynomial of order I and the 
coefficients a of the expansion (13) can be found as 


Ci := 


21+1 

2 

21+1 


R{tp) Pi {cos tp) sin tp dtp 


J \\ 


R{'tp)Pe{cos tp) sintp 


( d± 

\dX 


dX. 


(14) 


The integration limits Ai ,2 can be found from the condition tp{X) — 
0 and 'ip{X) — tv, respectively. Using this decomposition, it is 
straightforward to measure the differences between the value of the 
parametrized shadow at two different angles. For example, the rela¬ 
tive difference between the shadow atfi = 0 and at a generic angle 
Tp = TV/m can be computed simply as 


^ _R,p{Tp = 0) — R,p{Tp = Tv/m) 
Ri,{TP = 0) 

_ EfeoCrPr(cosV>)| ■j/>—7r/m 

Efco cePi{cosTp)\^=o 
When m = 1, expression (15) simplifies to 


5i := 1 - 


EEo(-l)^c^ 

E oo 
1=0 


(16) 


while, when m = 2 , the difference can still be computed analyti¬ 
cally and is given by 



(17) 


where we have introduced the following and more compact nota¬ 
tion that will be used extensively in the remainder 


^ = 0 ) = , (18) 
^■,=R^pf, = ^/2)=^yif,^^C2i, (19) 

OO 

^ ■.=R^{TP = ?,Tv/2) = Y^{-lYct. (20) 

1=0 


An areal radius is particularly useful as it enables one to 
measure two useful quantities, namely, the local deviation of the 
shadow R,p ~ R{Tp) from the areal circle, i.e., 

d,!, := \R^ - RA , ( 11 ) 

and its polar average 

1 1 

rf(y.> := ^ y d^dTp=—J \R^-RAdTp. (12) 

Note that although similar, the areal and the circumferential radii 
are in general different and coincide just for a spherically symmet¬ 
ric black hole, in which case R^ — Rc = R^, and of course 
d^ — 0 = ds- All of these geometrical quantities are shown 
schematically in Fig. 2. 


Analytic expressions for (15) when m > 2 are much harder to 
derive, but can be easily computed numerically. 

We note that the parametrization (15) is quite general and al¬ 
lows us to recover in a single definition some of the expressions 
characterizing the distortion of the shadow and that have been in¬ 
troduced by other authors. For example, the parameter can be 
associated to the distortion parameter S first introduced by Hioki & 
Maeda (2009) (cf. Fig. 3 of Hioki & Maeda 2009). Similarly, the 
parameter ^4 is directly related to the distortion parameter e intro¬ 
duced by Tsukamoto et al. (2014) (cf. Fig. 3 of Tsukamoto et al. 
2014). 

In what follows we will exploit the general expression for the 
polar curve representing the black hole shadow to suggest three 
different definitions that measure in a coordinate-independent man¬ 
ner the amount of distortion of the shadow relative to some simple 
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background curve, e.g., a circle. These expressions are all mathe¬ 
matically equivalent and the use of one over the other will depend 
on the specific properties of the observed shadow. 


parameter - (Ss.j can then be given by 


(Js.T : = 


ds, 

Us 


C2t- 


= : 


(29) 


3.1 Distortion parameter -1 

We start by considering three points on the polar curve A, B, and 
D, which occupy precise angular positions at t/i = 0, 7 r/ 2 , and 
37 r/ 2 , respectively (see diagram in Fig. 3). The corresponding dis¬ 
tances OA, OB and OD from the centre of coordinates O can then 
be expressed as 


= 0 ) = ^ ciPi{cosip)\-4,=o = , 


( 21 ) 


which reduces to the compact expression 

2^ 

Co 


(5s., 




(30) 


when only the first two coefficients in the expansion are taken into 
account, i.e., for co 7 ^ 0 7 ^ ci and ci = 0 , with £ ^ 2 (we recall 
that (5o,j = 0). We also note that the assumption (27) does not re¬ 
strict the analysis to spherically symmetric black hole space-times 
and, as we will show in Fig. 6 , the distortion parameter (29) can be 
applied also to axisymmetric space-times. 


Rg = 7r/2) = ^ C£Pr(cos V')lv,= 7 r /2 = SS, (22) 

1=0 

Rj^ :=R^{'<p = 3 tt/ 2) = Rg , (23) 

where in the last equality we have exploited the fact that the expan¬ 
sion in Legendre polynomials is symmetric with respect to the a 
axis. 

Next, we define a new parametric curve for which = 
R^ = R^ and thus that satisfies the following condition 


= .e/, 


or, equivalently, for £ > 0 


C2e-1 = C2i 


i-iY 


m 

22i(£\y 


- 1 


(24) 


(25) 


The corresponding polar expression, formulated in terms of 
the Legendre polynomials expansion (13), is therefore given by 


3.2 Distortion parameter - II 

A second possible definition of the distortion parameter is slightly 
more general and assumes that the radial distance of points A 
and B from the centre of coordinates is not necessarily the same, 
i.e., Rj^ y Rg. In this case, one can think of introducing a new 
point E on the a axis (this is shown schematically in the middle 
panel of Fig. Al, when considering the case of a Kerr black hole), 
such that the distances AE = EB and which could therefore serve 
as the centre of the reference circle. Since the values of the coor¬ 
dinates R^ and R^ are defined by expressions ( 21 ) and ( 22 ), we 
can use the condition AE = EB to find that one can easily find 
position of the point E on the a axis is given by 


Pr = 


pi 


2P. 


(31) 


with the corresponding angular position being either 0 or tt. 


Pv-,/ W = ^0 + X] 

t- 

P 2 l-l{cOS'ljj) + 


,J2£)^_ ■ 

' 22-*(£!)2 


P2r(cost/)) 


(26) 


To measure the distortion we need a reference curve, which 
we can choose to be the circle passing through the three points A, 
B, and D and thus with radius 


Ps.j ■- R^= SS = Rg = . 


(27) 


We can now compute the deviation of the parametric curve 
(26) from the corresponding background circle of radius at 
any angular position. However, as customary in this type of consid¬ 
erations, we can consider the shadow to be produced by a rotating 
black hole with spin axis along the /3 axis, so that the largest devia¬ 
tions will be on the axis of negative a (this is shown schematically 
in the left-hand panel of Fig. Al, when considering the case of a 
Kerr black hole). More specifically, we can define the difference 
between the curves at t/) = tt as 

00 00 

ds,i := Ps.i - R^,j {ip = tt) =SS “ ^ C2i-\ 


=2 C2r-i. 


(28) 


-i/lg = cos 


Pa - Pr 
|Pa - Pf 


(32) 


The radius of the circle passing through the three points A, B, and 
D is 


_ Rl+Rl 

2R, 


(33) 


so that the deviation of the shadow aX ip — it from the circle of 
radius Ps.jj can be found using the relation 


ds.ii — 2Ps,jj (P^ -F Re) ■ 


(34) 


Finally, we can introduce the distortion parameter ds,jj de¬ 
fined as 




(35) 


so that the second dimensionless distortion parameter is expressed 
as 




Po 


= 2 


.^ 2+^2 


(36) 


The expression for the dimensionless distortion (36) is in this 
case more complex that the one presented in Eq. (29); however, 
in the simpler case in which only the lowest order coefficients are 
retained, i.e., if cq 7 F 0 7 F ci and q = 0 for £ 2 , we have 


(5s.„ Ci 


2c? 


2cg -F 2coCi -F cf 


= 5i 


It follows that our first definition for the dimensionless distortion 


II • 


(37) 
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Figure 3. Schematic representation of the distortion parameter - III. The 
quantity measures the difference between the Legendre expanded 

poiar curve and the reference circie of radius ^nd passing 

through the points A, B, D, and having centre in point E. Aiso shown are 
the “zero-siope” points S and S'. 


3.3 Distortion parameter - III 

A third and possibly optimal definition of the distortion parameter 
is one that is meant to consider the case in which the shadow is still 
reflection symmetric relative to the a axis, but does not cross the /3 
axis with a zero slope. Rather, the curve admits a point, say S, at 
angular position 0 < < tt, where it has zero slope relative to 

the (a, j3) coordinate system (see diagram in Fig. 3 and the right- 
hand panel of Fig. A1 for the case of a Kerr black hole). This point 
will be referred to as the “slope point” of the parametric curve 
representing the shadow. 

To compute the position of this point in the {a, j3) coordinates 
we simply need to find the solution to the equation 


dot 


'I's 


= 0 , 


or, equivalently, solve for the differential equation 


dRip 

dtp 


sin tp + Rp; cos tp = 0 . 


(38) 


(39) 


Using the expansion in terms of Legendre polynomials (13), 
we can rewrite equation (39) as 


cePi{x)x - (1 - = 0, (40) 


dx 


where we have set x := costp. The solutions of (40) provide the 
positions of all the possible slope points in the parametric curve, 
and the solution is unique in the case in which the shadow R{tp) is 
convex. The corresponding coordinates of the point S are then 


OO 

Rs =y2cePt{Xs), (41) 

£^0 

tp^ =cos~^{Xg). (42) 

As for the second distortion parameter in Sect. 3.3, we set E 
to be the centre of the circle passing through the points A, S, and 


S', where S' the point is symmetric to the point S with respect to 
the a axis. Using the condition AE = ED, we obtain the solution 


Rr, = 


tp^ = cos 


2 EfcoC<?(l - Pdxs)^s) 

R, 


\Rr 


(43) 


(44) 


Also for this third case, the radius of the circle Rs^m pass¬ 
ing through the three points A, S, and S', the distortion parameter 
ds.jjj, and its dimensionless version ds.m, have respectively the 
form 


Rs.m — 


‘iT.ZoCiil - PdXs)Xs) 

'^s.iii ~ — {Rji A Rc) ~ I y ] cePi{Xg) 

\r=o 

(E^o cePt{Xs) - Xs E^i C2r-i - V) 

Y.Z,cdl-PdXs)Xs) 

1 / OO 

'5^.111 = = 2 ij^ciPiixA 

(Efco DPi{Xs) - Xs Efci C 2 r -1 - 


(45) 


(46) 


(47) 


_ 2XsJ^J2T=oDPeiXs) -h iJ2T=o^tPdXs)d ' 

We note that this definition is similar to the one proposed by 
Hioki & Maeda (2009), who measure the dimensionless distortion 
of the shadow as 

ds 


X _ 


(48) 


where 


d-s,am ■= Rsid = tt) - Rs,nm y ( 49 ) 

and with i?s,HM being the radius of the circle passing through the 
points A, S, and S'. The most important difference with respect 
to the definition of Hioki & Maeda (2009) is that we here express 
the parametric curve in terms of the general Legendre expansion 
(13), while Hioki & Maeda (2009) assume the knowledge of Rp, at 

tp = TT. 

Also in this case, expressions (45)-(47) are not easy to handle 
analytically. However, in the simplest case in which the expansion 
(13) has only two nonvanishing terms, such that cq 7 ^ 0 7 ^ ci and 
Cl = 0 for £ ^ 2, Eq. (40) takes the more compact form 

2cia:^ -|- cox — ci = 0 , (50) 

with solution 



(51) 


and where the -|- or — signs refer to when ci > 0 and ci < 0 , 
respectively. The corresponding quantities Rs.^i, ds,ni (5s,jj; 

have in this case the following form 


_2co “t“Ci(l 4“ 2coCi (1 A x^') 

2[co 4“ Ci(l 4“ 

_ Ci(l4-Xg) 

Co + Ci(l -j- Xg) 

_ _ 2ci(l 4- Xg) _ 

2co 4-Ci(l 4-Xg) 4-2coCi(l 4-Xg) 


(52) 

(53) 

(54) 
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If the shadow is perfectly circular with radius cq, then ci = 0 
and expressions (52)-(54) show that -Rs.m = cq, = 0 = 

as expected. 


4 APPLICATION OF THE FORMALISM TO BLACK 
HOLE SPACE-TIMES 

Having constructed a general formalism that allows us to describe 
in a coordinate-independent manner the black hole shadow and to 
measure its deformation, we are now ready to apply such a formal¬ 
ism to the specific case of some well-known space-time metrics 
referring to axisymmetric black holes. In particular, we will obvi¬ 
ously start with the application of the formalism to a rotating (Kerr) 
black hole (in Sect. 4.1), to move over to a Bardeen black hole and 
to a Kerr-Taub-NUT black hole in Sect. 4.2. We note that we do 
not consider these last two examples of black holes because they 
are particularly realistic, but simply because they offer analytic line 
elements on which our formalism can be applied. 


4.1 Kerr black hole 


We start with the Kerr space-time, whose line element in Boyer- 
Lindquist coordinates reads 


ds = - 1 - 


2Mr 




^ dC - 


4aMr sim 6 


dtd(f> 


, 2 2 2a‘‘Mrshr6 

+ \ C + a^ + ^ 


+ ^dr^ + T.de'^ . 

A 


sin^ 9d(f 


where 


(55) 


_ 2 . 2 2/3 

:= r -I- a cos 0 , 


A — 2Mr + . 


(56) 


with M being the mass of the black hole and a := J/M its specific 
angular momentum. 

Since the shape of the shadow is ultimately determined by 
the innermost unstable orbits of photons, hereafter we will con¬ 
centrate on their equations for photons. In such a space-time, the 
corresponding geodesic equations take the form 


dt'\ AE — 2aMrLz 

^ A 



_ 2aAdrE; (S - 2Mr) Lz 
\ dA / A A sin^ 0 


(57) 

(58) 

(59) 

(60) 


where A is an affine parameter, 

7^ ■-E'^r* + {a^E^ - LI - Q) 

+ 2M \{aE — Lzf -I- Q] r — a^Q , (61) 

0 —Q (a^E^ — L\ csc^ 9) cos^ 9 , (62) 

A := — a^Asin^ 9 , (63) 

and E and Lz are the photon’s energy and the angular momentum, 
respectively. The quantity Q 


Q=Pb + cos^ 9 



2 


— a 



(64) 


is the so-called Carter constant and pg ;= T,d9/dX is the canonical 
momentum conjugate to 9. 

Using these definitions, it is possible to determine the unstable 
orbits as those satisfying the conditions 


7^(f) 


dTZjf) 

dr 


and 


d^TZjf) 

dr 2 


>0, 


(65) 


where f is the radial coordinate of the unstable orbit. Introducing 
now the new parameters ^ := LzjE and 77 := Q/E'^, the celestial 
coordinates a and /3 of the image plane of the distant observer are 
given by (Bardeen et al. 1972) 


a = , P = ±(77 -I- o' cos^ i — cot^ , ( 66 ) 

sm 7 

where i is the inclination angle of the observer’s plane, that is, the 
angle between the normal to the observer’s plane and the black 
hole’s axis of rotation. In the case of the Kerr space-time (55) and 
after using the conditions (65), one can easily find that the values 
of ^ and 77 relative to the circular orbit (c) are (Young 1976; Chan¬ 
drasekhar 1998) 


= 


M{r^ — a^) — r(r^ — 2Mr + a^) 
a(f — M) 

a’^{f — M)2 ’ 


(67) 

( 68 ) 


Next, to investigate the shape of the black hole shadow we 
introduce the generic celestial polar coordinates (7?', 7 /)') (cf. Sect. 
2 and Fig. 1) defined as 


R! ={a , 

(69) 


(70) 


Assuming for simplicity that the observer is in the equatorial plane, 
i.e., that i — 7r/2, then in terms of the {R'jtp') coordinates the 
shadow of black hole can be described as (hereafter we will set 
M = 1) 


,, (2r* + 2a?r — 6r^ + o? + aPr'^f^^ 

r — 1 

/ ^ -1 ( {r^[4a^-r(r-3)^]}^/^ 

\r2 — a2 — r{r2 — 2r + a^) 


(71) 

(72) 


Making use of the procedure described in Sect. 2, it is straight¬ 
forward to determine the coordinates (3)-(4) of the effective centre 
of the shadow, and to perform the Legendre expansion (13). To the 
best of our knowledge, no analytic expression exists to cast the co¬ 
ordinates (71)-(72) as a polar curve = R{tp)- However, such a 
curve can be easily constructed numerically and from it the Legen¬ 
dre expansion (13) can be computed. 

Figure 4 summarizes the results of our approach by report¬ 
ing in the left-hand panel and in a logarithmic scale the normalized 
values of the expansion coefficients ci as a function of the Legen¬ 
dre order £. Different curves refer to the different values consid¬ 
ered for the dimensionless spin parameter a* := J/M'^ = a/M, 
which ranges from a* = 0.4 (blue solid line) to a* = 0.99 (red 
solid line). Interestingly, the series converges very rapidly (essen¬ 
tially exponentially) and already with f = 4, the contribution of 
higher-order terms is of the order of 10“^, decreasing further to 
~ 10“^ for 1 = 6. Furthermore, even when considering the more 
severe test of a* = 0.99, the expansion coefficient with 1 = 6 
is only a factor 2-3 larger than the corresponding coefficient for a 
slowly rotating black hole. The right-hand panel of Fig. 4 shows a 
direct measure of the relative differences between the polar curve 
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Figure 4. Left-hand panel: Magnitude of the expansion coefficients of the polar curve and shown as a function of the expansion order. Note the very 
rapid (exponential) convergence of the expansion. The coefficients are computed for a Ken' black hole and different lines refer to different values of the spin 
parameter a*. Right-hand panel: Relative differences between the polar curve for the black hole shadow as constructed from expressions (71)-(72) and the 
corresponding curve obtained from the expansion, i.e., 1 — ciPijR^. Different lines refer to different truncations of the expansion and show that three 
coefficients are sufficient to obtain deviations of a few percent. 





0123456789 10 
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Figure 5. Magnitude as a function of the expansion order of the three different distortion parameters defined as in Eqs. (29), (36), and (47), thus 

measuring the difference between the Legendre expanded polar curves and the reference circles of radii All curves refer to a Kerr 

black hole and different colours are used to represent different values of the spin parameter. 


for the black hole shadow as constructed from expressions (71)- 
(72) and the corresponding curve obtained from the expansion, 
i.e., 1 — ciPi/R^. Remarkably, already when considering the 
first three terms in the expansion, i.e., co,Ci, and C2, the relative 
difference is of a few percent (blue line), and this further reduces to 
10“^ when the expansion is truncated at n = 4 (black line). 

In summary. Fig. 4 demonstrates that when considering a Kerr 
black hole, the approach proposed here provides a coordinate inde¬ 
pendent and accurate representation of the black hole shadow and 
that a handful of coefficients is sufficient for most practical pur¬ 
poses. In the following Sections we will show that this is the case 
also for other axisymmetric black holes. 

Before doing that, we show in Fig. 5 the values of the dimen¬ 
sionless distortion parameters as computed for the shadow of a Kerr 
black hole and for increasing values of the expansion index 1. The 
three different panels are relative respectively to the parameters 


(29), (36), and (47), with the different curves referring to values 
of the spin parameter a*, ranging from a* = 0.4 (blue solid line) 
to a* = 0.99 (red solid line). As one would expect, for all values of 
a*, each of the three distortion parameters decreases as the expan¬ 
sion includes higher-order terms. At the same time, because larger 
rotation rates introduce larger distortions in the shadow, they also 
lead to larger values of the distortion parameters for a fixed value 
oft 

Finally, Fig. 6 offers a comparative view of the different dis¬ 
tortion parameters for specific values of the spin parameter, with 
the left-hand and right-hand panels referring to a* = 0.4 and 
a* — 0.99, respectively. This view is rather instructive as it shows 
that the different definitions lead to significantly different values 
of the distortion, despite they all refer to the same parametric polar 
curve. Furthermore, it helps appreciate that the distortion parameter 
5s,jj is systematically smaller than the other two and hence not the 
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Figure 6. Comparative view of the different distortion parameters (red line), 5s,jj (blue line), and 5s,jjj (light-blue line). The left-hand and right-hand 
panels show the values of the distortion parameters as a function of the expansion order (cf. Fig. 5), and refer to a Kerr black hole with spin a* = 0.40 and 
a* = 0.99, respectively. 


optimal one. This is because a larger value of the distortion param¬ 
eter will increase the possibility of capturing the complex structure 
of the shadow. The fact that the curves for 5s,i and 5s,intersect 
for the Kerr black hole considered at £ = 5 implies that both dis¬ 
tortion parameters (29) and (47) are useful indicators, with 5s,m 
being the recommended choice for expansions having £ ^ 5. 

4.2 Bardeen and Kerr-Taub-NUT black holes 

We continue our application of the formalism developed in Sects. 2 
and 3 by considering the space-time of a rotating Bardeen black 
hole (Bardeen 1968). We recall that in Boyer-Lindquist coordi¬ 
nates, the metric of a Kerr and of a Bardeen black hole differ 
uniquely in the mass, which needs to be modified as (Bambi & 
Modesto 2013; Tsukamoto et al. 2014) 


hole with nonvanishing gravitomagnetic charge n and specific an¬ 
gular momentum a := J/M. The corresponding metric is given 
by (Newman et al. 1963) 

ds^ = — i (A — sin^ 0) dC + E + dO^'^ 

+ i [(E + axf sin^ 6 - x^A] d(f 
+ —a.{T,-\-axitAtd 6)dtd4>, (78) 

where the functions A, E, and x are now defined as 
A :=r^ + — rd — 2Mr , 

E :=r^ + (n + a cos Of , 

X :=a sin^ 0 — 2ncos 0 . (79) 


/ ^2 y/2 

= , (73) 


where the parameter g is the magnetic charge of the nonlinear elec¬ 
trodynamic field responsible for the deviation away from the Kerr 
space-time. 

The impact parameters ^ and g relative to the circular orbit are 
in this case (Tsukamoto et al. 2014) 


m[(2 — f)r^ — fa^] — r{r^ — 2mr + a?) 
a{f- fm) 

r^{4(2 — f)a^m — f[r — (4 — f)m]^} 
a^{f — fmf ’ 


(74) 

(75) 


and can be taken to define the shadow of black hole. Note that m 
and / are functions of the unstable circular radius f and are given 
by 


m =m(r) = M 




3/2 




(76) 


(77) 


j,2 _|_ g2 

In complete analogy, we can consider a Kerr-Taub-NUT black 


In this case, the impact parameters ^ and rj for the circular 
orbits are given by (Abdujabbarov et al. 2013) 


^ _a^(l + r) + r^(r — 3) + rd{l — 3r) 

~ a(l - f) ’ 

7)c = ^ —— |f®[4a^ - f(r - 3)^] - rd [4r^a^ 

a‘‘{r — 1 ) [ 

+ (1 - 3T)(n^(l - 3r) - 6f^ + 4aV + 2f®)] I , 


(80) 


(81) 


and define the shadow of the Kerr-Taub-NUT black hole. 

Applying the formalism described in Sect. 2, it is possible to 
compute the coefficients of the Legendre expansion (13) also for a 
Bardeen and for a Kerr-Taub-NUT black hole. The numerical val¬ 
ues of these coefficients are reported as a function of the expansion 
order £ in Fig. 7, where the left-hand panel refers to a Bardeen 
black hole, while the right-hand one to a Kerr-Taub-NUT black 
hole. More specifically, the different lines in the left-hand panel re¬ 
fer to different values of the magnetic charge: g = 0.00 (red line), 
g — 0.30 (blue line), and g = 0.60 (light-blue line); all lines refer 
to a fixed value of the rotation parameter a* = 0.60. Very similar is 
also the content of the right-hand panel of Fig. 7, which is however 
relative to a Kerr-Taub-NUT black hole. The different curves in this 
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Figure 7. Left-hand panel: Magnitude of the expansion coefficients as a function of the expansion order I for the different values of the magnetic charge 
of a Bardeen black hole: g = 0.00 (red line), g = 0.30 (blue line), and g = 0.50 (light-blue line). All lines refer to a fixed value of the rotation parameter 
a* =0.60 (cf. left panel of Fig. 4). Right-hand panel: The same as in the left-hand panel but for a Kerr-Taub-NUT black hole. Different curves refer to the 
different values of the NUT parameter: n = 0.00 (red line), n = 0.30 (blue line), and n = 0.60 (light-blue line). All lines refer to a fixed value of the rotation 
parameter a* = 0.60. 


case refer to the different values of the NUT parameter: n = 0.00 
(red line), n = 0.30 (blue line), and n = 0.60 (light-blue line); 
once again, all lines refer to a fixed value of the rotation parameter 
a* = 0.60. 

In analogy with what shown in the left-hand panel of Fig. 4 
for a Kerr black hole, also for these black holes the series con¬ 
verges very rapidly and already with 1 = 4, the contribution of 
higher-order terms is of the order of 10“^, decreasing further to 
~ 10~® for £ = 6, even when considering higher larger mag¬ 
netic charges or NUT parameters. Furthermore, in analogy with the 
right-hand panel of Fig. 4, we have checked that the relative differ¬ 
ences between the polar curves for the shadow constructed from 
expressions (74)-(75) and (80)-(81), and the corresponding curve 
obtained from the expansion, i.e., 1 — CiPi/R- 4 >, is below 10“^ 
when n = 2 (not shown in Fig. 7); this difference further reduces 
to 10“"* when the expansion is truncated at n = 4. 

In conclusion, also Fig. 7 demonstrates that the approach pro¬ 
posed here provides a coordinate independent and accurate repre¬ 
sentation of black hole shadows in space-times other than the Kerr 
one. 


5 COMPARISON WITH NOISY OBSERVATIONAL DATA 

All of the considerations made so far have relied on the assumption 
that the shadow is a well-defined one-dimensional curve (cf. dis¬ 
cussion in Sect 2). In practice, however, this is not going to be the 
case as the astronomical observations will have intrinsic uncertain¬ 
ties that, at least initially, will be rather large. It is therefore natural 
to ask how the formalism presented here will cope with such un¬ 
certainties. More precisely, it is natural to ask whether it will still 
be possible to determine the effective centre of a noisy polar curve 
and then determine from there its properties. Although a very ob¬ 
vious and realistic problem, this concern is systematically ignored 
in the literature, where the shadow is traditionally assumed to have 
no uncertainty due to the observational measurements. 

While awaiting for actual observational data, we can straight¬ 


forwardly address this issue in our formalism and mimic the nois¬ 
iness in the observational data by considering the polar curve as 
given by the Legendre expansion (13), where however the differ¬ 
ent coefficients ci are artificially perturbed. More specifically, we 
express the shadow via the new expansion 

OO 

+ Ac)P^(cosV^), (82) 

^-0 

where Ac is a random real number chosen uniformly in the range 
[—Amax, Amax]. In tMs Way, our putative polar curve represent¬ 
ing the shadow will be distorted following a random distribu¬ 
tion and we have optimistically assumed a variance of only 5%, 
i.e., Amax = 0.05. Of course, there is no reason to expect that the 
error distribution in the actual observational data will be uniform, 
but assuming a white noise is for us the simplest and less arbitrary 
choice. 

With the setup described above and the formalism discussed 
in the previous Sections, we have considered a reference shadow of 
a Kerr black hole with spin parameter a/M = 0.99 and have repro¬ 
duced it after truncating the expansion (82) at £ = 9, which is more 
than sufficient given the accuracy obtained at this order (cf. Fig. 4). 
We have therefore constructed a very large number of such realiza¬ 
tions of the observational shadow after making use of Atot = 10® 
draws of the random deformation Ac. For each putative observed 
reconstructed shadow we have computed the distortion parameters 
5s,I — defined in Eqs. (29), (36), and (47), as well as the dis¬ 
tortion definition of Hioki & Maeda (2009) and defined in Eq. (48). 

Eor each of the shadow realizations we have therefore com¬ 
puted the measurement error as 

e* := 5s - 5s.* , (83) 

where 5s is the exact distortion of the background Kerr solution 
and measuring the relative difference of the shadow at t/; = 0 and 
t/) = TT. On the other hand, 5s,* is given by either 5s, j — 5s, or 

5s,HM ■ 

Figure 8 shows the distributions of the errors computed in this 
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Figure 8. Comparison of probability density distributions of the errors e* 
computed in the measurement of the distortion parameters 5^ j — <5s,ni 
a Ken' black hole shadow reconstructed using a perturbed expansion [cf. Eq. 
(82)]. Also shown is the distribution of the error measured when using the 
distortion parameter introduced by Hioki & Maeda (2009) and that has a 
larger variance. 


way for the four different possible definitions of the distortion pa¬ 
rameters, with the black line referring to the distortion parameter 
in Eq. (48), and the red, blue and light-blue lines referring to the 
definitions (29), (36), and (47), respectively. Note that the values 
of the probability densities distributions are reported in such a way 
that 



The distributions reported in Fig. 8 are rather self-explanatory. 
Clearly, all the different definitions are centred on e* = 0, indicat¬ 
ing that on average they provide a good measurement of the distor¬ 
tion. On the other hand, the variance of the different distribution is 
rather different. Overall, the distortion parameters — (5s,m have 
comparable variances, with a slightly smaller variations for the def¬ 
inition 5s,i - However, the variance of the distortion parameter for 
5s,HM almost twice as large as the others and it essentially spans 
the 5% variation that we introduce in the random distortion Ac. 
These results are rather reassuring as they indicate that new defini¬ 
tions are not only accurate, but also robust with respect to random 
white noise. Furthermore, they appear to be superior to other dis¬ 
tortion measurements suggested in the past. 

As a final remark we note that the introduction of the pertur¬ 
bations in the expansion (82) also has the effect of changing the 
position of the effective centre of the shadow and hence the values 
of Rq and i/jq [cf. Eqs. (2) and (4)]. Fortunately, such variations 
represent only a high-order error, which is much smaller than those 
measured by the distortion parameters, with maximum measured 
variance of the order of 10“^. As a result, the distortions reported 
in Fig. 8 are genuine measurements of the shadow and not artefacts 
introduced by the changes in the effective centres. 


6 CONCLUSION 

The radio-astronomical observations of the shadow of a black hole 
would provide convincing evidence about the existence of black 
holes. Further, the study of the shadow could be used to learn about 
extreme gravity near the event horizon and determine the correct 
theory of gravity in a regime that has not been explored directly so 
far. A number of different mathematical descriptions of the shadow 
have already been proposed, but all make use of a number of sim¬ 
plifying assumptions that are unlikely to be offered by the real ob¬ 
servational data, e.g., the ability of knowing with precision the lo¬ 
cation of the centre of the shadow. 

To remove these assumptions we have developed a new gen¬ 
eral and coordinate-independent formalism in which the shadow is 
described as an arbitrary polar curve expressed in terms of a Leg¬ 
endre expansion. Our formalism does not presume any knowledge 
of the properties of the shadow and offers a number of routes to 
characterize the properties of the curve. Because our definition of 
the shadow is straightforward and unambiguous, it can be imple¬ 
mented by different teams analysing the same noisy data. 

The Legendre expansion used in our approach converges ex¬ 
ponentially rapidly and we have verified that a number of coeffi¬ 
cients less than ten is sufficient to reproduce the shadow with a 
precision of one part in 10®, both in the case of a Kerr black hole 
with spin parameter of a/M = 0.99, and in the case of Bardeen 
and Kerr-Taub-NUT black holes with large magnetic charges and 
NUT parameters. Furthermore, the use of a simple Legendre ex¬ 
pansion has allowed us to introduce three different definitions of 
the distortion of the shadow relative to some reference circles. The 
comparison of the different definitions has allowed us to determine 
which of them is best suited to capture the complex structure of the 
shadow and its distortions. 

Finally, again exploiting the advantages of the Legendre ex¬ 
pansion, we were able to simulate rather simply the presence of 
observational random errors in the measurements of the shadow. 
Constracting a large number of synthetically perturbed shadows, 
we have compared the abilities of the different parameters to mea¬ 
sure the distortion in the more realistic case of a noisy shadow. 
Overall, we have found that our new definitions have error distri¬ 
butions with comparable variances, but also that these are about a 
factor of 2 smaller than the corresponding variance measured when 
using more traditional definitions of the distortion. Given these re¬ 
sults, the approach proposed here could be used in all those studies 
characterizing the shadows of black holes as well as in the analysis 
of the data from experimental efforts such as EHT and BHC. 
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APPENDIX A: GRAPHICAL REPRESENTATION OF 
DISTORTION PARAMETERS: KERR BLACK HOLES 

To help in the visualization of the different distortion parameters 
introduced in Section 3, we here present their graphical represen¬ 
tation when they are applied to the shadow of a Kerr black hole 
shadow with spin parameter a* = 0.98. The left-hand panel of 
Fig. A1 represents the distortion parameter I in Sect. 3, where the 
reference circle is defined in a such a way that the points B, D 
and the centre of the reference circle are all on the vertical axis j3. 
Hence, the difference of the radial coordinates of the point C of the 
shadow boundary and of the left-hand point of the circle intersect¬ 
ing the axis a corresponds to the distortion parameter dsj. 

Similarly, the middle panel of Fig. A1 refers to the distortion 
parameter II, where the reference circle passes through the points 
A and B, which are on the axis /J. The centre of the reference cir¬ 
cle is on the point E and does not coincide with the centre of the 
coordinate system. The position of the reference circle centre E is 
instead defined by Eq. (31). The difference of the radial coordinates 
of the point C of the shadow and of the left-hand point of the circle 
intersecting with the axis a corresponds to the distortion parameter 
ds,ii. 

Finally, the right-hand panel of Fig. A1 corresponds to the dis¬ 
tortion parameter III, where the reference circle passes through the 
slope points S, S', and the right-hand point A of the shadow in¬ 
tersecting the axis a. The positions of the slope points are defined 
by solving Eq. (39). The centre of the reference circle does not co¬ 
incide with the origin of the coordinate system and its position is 
defined by Eq. (43). The difference of the radial coordinates of the 
point C of the shadow and of the left-hand point of the circle inter¬ 
secting the axis a corresponds to the distortion parameter ds,iii. 

Note that the radius of the reference circle radius depends on 
the definition used for the distortion parameter; hence the distortion 
parameters are also different for each definition. 
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Figure Al. Schematic representations of the distortion parameters 1,11 and III when applied to the shadow of a Kerr black hole. The left-hand panel refers to 
the definition I, where the centre of the reference circle as well as the points B and D (which are not slope points) are on the coordinate axis The middle 
panel refers instead to definition II and in this case the centre of the reference circle E is displaced along the a axis. Finally, the right-hand panel shows the 
representation of the definition III, where we consider the reference circle passing through the point A and the slope points B and D are not on the (3 axis. 
The centre of the reference circle E is also displaced and not exactly at the centre of the coordinate system. 












